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Abstract 

We investigate constraints on the time variation of the fine structure constant between the 
recombination epoch and the present epoch, Aa/a = (a rec — «now)/anow, from cosmic microwave 
background (CMB) taking into account simultaneous variation of other physical constants, namely 
the electron mass m e and the proton mass m p . In other words, we consider the variation of 
Yukawa coupling and the QCD scale Aqcd in addition to the electromagnetic coupling. We clarify 
which parameters can be determined from CMB temperature anisotropy in terms of singular value 
decomposition. Assuming a relation among variations of coupling constants governed by a single 
scalar field (the dilaton), the 95% confidence level (C.L.) constraint on Aa/a is found to be 
—8.28 x 10~ 3 < Aa/a < 1.81 x 10~ 3 , which is tighter than the one obtained by considering 
only the change of a and m e . We also obtain the constraint on the time variation of the proton- 
to-electron mass ratio fi = m p /m e to be —0.52 < A^/u < 0.17 (95% C.L.) under the same 
assumption. Finally, we also implement a forecast for constraints from the PLANCK survey. 
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I. INTRODUCTION 



Whether or not the physical constants are truly constant is not only a fundamental is- 
sue in physics but also an important probe of the theories of extra dimensions. Since the 
early studies of Dirac , many theoretical models that could accomodate the variation 
of physical constants have been proposed . Unification theories such as super- 

string theories predict the existence of an additional scalar field $ called "dilaton" and the 
destabilization of this field j^l, 9] might cause the variation of coupling constants. In more 
phenomenological contexts, dynamical scalar fields introduced to explain the recent discovery 



fields in the S 
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of the accerelated expansion may be non-minimally coupled to gauge 
model, which may also lead to time varying coupling constants [lOl . 
Thus, investigating the variation of the coupling constants can be an important probe of 
these theories. 

Not only these theoretical sides but also the observational analysis of various systems 
have been playing an important role in the investigation of the time variation of physical 
constants. Among them, the claim for the deviation of the past value of the fine structure 
constant a from the present one by using the high-redshift quasar absorption system 
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had the great impact on this topic. Further analysis for the quasar absorption s 



the support for this positive result of a cosmologica 
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variation of a 
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2l|, though 
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some authors takes objection to such variations 

Stimulated by these studies from both theoretical and observational sides, various con- 
straints on the time variation of the fine structure constant a have been investigated from 
diverse observations. We briefly summarize those terrestrial and celestial limits on a as 
follows (see for a review). The atomic clocks constrain the current value of the 
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. The measure- 



temporal derivative of a as a/a = (—3.3 ± 3.0) x 10~ 16 yr~ x 25|, |26|, 
ment of the frequency ratio of aluminium and mercury single-ion optical clocks provides 
a/a = (—1.6 ± 2.3) x 10~ 17 yr _1 {2^. From the analysis of Sm isotopes in the Oklo natural 
reactor in Gabon, we get two bounds on the variation of a as Aa/a = —(0.8 ± 1.0) x 10 -8 
and Aa/a = (0.88 ± 0.07) x 10~ 7 [30|, which measures the value at the redshift z ~ 0.1 



constraint. From the spectra of quasars, several limits have been obtained using vari- 



ous data sets: Aa/a = 
instrument [18], Aa/a 



'-0.57 ± 0.11) x 10" 



[z ~ 



-0.64 ± 0.36) x 10~ 5 



0.2 - 4.2) from the Keck/HIRES 
~ 0.4 — 2.3) from the Ultravio- 
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let and Visual Echelle Spectrograph (UVES) instrument [19|, l20j. Big bang nucleosyn- 
thesis (BBN) provides constraints at very hi gh redshifts (z ~ 10 9 — 10 10 ), for example, 
-5.0 x 1(T 2 < Aa/a < 1.0 x 10~ 2 (95%C.L.) [31|. Finally Cosmic Microwave Background 
(CMB) measures a at z ~ 10 3 and the limits from WMAP 1-year, 3-year and 5-year data 
read -0.06 < Aa/a < 0.01 [2a], -0.039 < Aa/a < 0.010 |34j and -0.028 < Aa/a < 0.026 



351 ] respectively, all of which are at 95%C.L. Recently, Menegoni et al. [36] obtained more 



improved result —0.013 < Aa/a < 0.015 at 95%C.L. from WMAP 5-year date combined 
with ACBAR, QUAD and BICEP experiments data. 

Along with a, several observational constraints have been obtained on another physical 
constant, the proton-to-electron mass ratio /i = m p /m e . As in a, atomic clocks very tig htly 
constrain the current value of its temporal derivative as fi/fj, — (— 1.5± 1.7) x 10~ 15 yr _1 25]. 
Reinhold et al. [3^] reports the non-vanishing variation of /i as A/i//i = (2.4 ± 0.6) x 10~ 5 
from a weighted sum of accurate H2 spectral lines in two quasar systems (corresponding 
redshifts are z = 2.59 and z = 3.02). Based on the quasar absorption spectra of NH 3 , the 
limit is Afi/fi = (0.6 ± 1.9) x 10~ 6 at z = 0.6847 [38;] and further detailed measurement by 



the same method tighten the bound as A/i//i = (0.74 ± 0.47) x 10~ 6 39]. 

In this paper, using WMAP 5-year data, we focus on the CMB constraint on a in the 
case that the multiple physical constants may vary with time simultaneously according to 



the expectation values o 
for BBN 
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a single dilaton field Although similar analyses has been done 
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44j , quasar absorption system 45j, |46j or the Oklo reactor 



47], few attempts from CMB data have been made except for several preceding literatures 
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49J, which treated simultaneous variations of only a and m e . In addition to a 



and m e , many theoretical models predict the time variation of the QCD scale Aqcd (which 
results from the variation of the strong coupling constant), so it is important to include 
the Aqcd effect in the analysis. We consider the effect of varying Aqcd as the change in 
the proton mass m p which is assumed to be proportional to Aqcd 1 • Thus, we consider 
simultaneous variation of a, m e and m p . Since how the variation of Aqcd is related to that 
of a or m e depends on the details of the unification model [50], we here adopt a specific 



1 When we consider the time variation of Aqcdj the neutron mass also changes. In this paper, we neglect 
the difference between the proton mass and the neutron mass, and when we refer to the time variation 
of the proton mass m pi we implicitly include the time variation of the neutron mass assuming that the 
neutron mass changes in exactly the same way as the proton mass. 
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example of the relation among those physical constants following 31] or [40|], which is based 
on a low energy effective theory of string theory. The newly incorporated change of Aqcd 
leads to a new and tighter constraint on the time variation of a. We also obtain a constraint 
on /i induced by the time dependence of Aqcd and m e . 

This paper is organized as follows. Section [III provides how the variation of each phys- 
ical constant affects the CMB power spectrum both qualitatively and quantitatively with 
particular attention to Aqcd- Section [TTT1 briefly describes the model we adopt based on 
string theory which governs the time evolution of the three physical constants we consider, 
namely, the variation of the fine structure constant, the electron mass and the proton mass. 
In section I FY] we provide the limits of those variations. Section |V] is devoted to conclusion. 
The details of the calculations of the degeneracy and parameter dependence are given in 
the appendix A. In appendix B, we present a forecast for constraints from the PLANCK 
experiment. 

II. THE EFFECT OF VARIATIONS OF PHYSICAL CONSTANTS ON THE CMB 

As is well known, a primary CMB photon is a relic signal from the last scattering surfuce 
about thirteen billion years ago, which means that the most important process for CMB 
photon is the recombination. It has been known that changing the value of the fine structure 
constant a and the electron mass m e affects the CMB power spectrum mainly through the 
change of the epoch of recombination 
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54] . The larger value of a or m e at 



the recombination epoch causes the higher redshift of the last scattering surface because the 
hydrogen binding energy scales as a 2 m e . This results in three characteristic signatures in the 
angular power spectrum of the temperature anisotropy, namely, shift of the peaks to higher 
multipoles, increase of the height of the peaks due to the enhanced early integrated Sachs 
Wolfe effect, and decrease of the small-scale diffusion damping effect. To be more precise, 
however, the process of the recombination depend on a and m e in a subtly different manner 
as will be explicitly shown in eqs.Q-Q in the argument of the four CMB-characteristic 
variables below. 

The time variation of the proton mass m p affects the CMB power spectrum mainly 
through the baryon density p#. In the standard cosmology with no variation of physical 
constants, the baryon density scales as pB,rec = AB,o(«o/«rec) 3 , where a is the scale factor and 
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the subscripts "rec" and "0" denote values at the epoch of recombination and the present 
epoch, respectively. When we allow the time variation of m p , there is an additional change 
as 

\ 3 / \ 3 / 

ao \ m P ,rec __ f a \ ( 1 , Am, 



pB,rec = PB,0 = PB,0 I H , (1) 

\a rcc J m pfi V a rec/ V m v J 

where we define Am p = m p ^ ec — m Pj o- Therefore, the effect of changing m PiI . ec is very 
similar to changing p B $. One may even consider that m PyTec and p B ,o are totally degenerate. 
However, there appears a subtle difference as regards the baryon number density % which 
affects the recombination process. This is because, contrary to the case of p B , the scaling 
^s.rec = ^B,o( a o/ a rec) 3 is not affected by the variation of m p (note that n B = p B /m p ). Thus, 
suppose one increases p B<0 an d decreases m PjTec while preserving the value of pB,rec, ^B,rec is 
increased. This is the reason why varying p#,o (as is done in the usual parameter search) and 
varying m p ^ ec are different and we can break the degeneracy between them in the parameter 
estimation although the degeneracy turns out to be rather strong. 

With the careful attention to the difference between p B and Am p noted in the last para- 
graph, we calculated the angular power spectra of the CMB temperature anisotropy, Cf T , 



for several different values of Am p using the modified CAMB code [55|]. The result is shown 
in FIG. [U where we also draw the cross correlation of temperature and E-mode polariza- 
tion, Cf E , and the angular power spectrum of E-mode polarization, Cf E , for completeness. 
Incorporating Am p into the numerical calculation, we assume m p to change instantaneously 
at a specific time (z = 100). In other words, m p ^ ec was constant before z = 100 and after 
that the proton mass becomes the present value 2 . This procedure is adopted throughout 
this paper. We consider a flat ACDM universe with a power-law adiabatic primordial fluc- 
tuation and take the fiducial standard cosmological parameter values as the WMAP 5-year 
marginalized mean values 
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{VL B h 2 , VL DM h 2 , h, r, n s ) = (0.0227, 0.1099, 0.72, 0.087, 0.963), (2) 

where Q B is the baryon energy density p B $ normalized by the critical density p C rit,cb ^dm is 
the dark matter energy density Pdm,o normalized by the critical density p crit , h is the current 



2 Here, only the mass difference Am p matters and it is not important how the proton mass has changed 
between the recombination epoch and the present epoch, because characteristic signatures such as acoustic 
oscillation or Silk damping are imprinted before and during the recombination process, which is also true 
on the effects that Am. n induces. 
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FIG. 1: The top left panel is the CMB temperature anisotropy spectra for no change of m p (solid 
red curve), an increase of m p by 10% (dashed green curve), a decrease of m p by 10% (dotted blue 
curve). The middle left and bottom left panels are the TE and EE spectra respectively for three 
different values of Am p /m p . For comparison, the corresponding spectra for three different values 
of Qsh 2 are shown in the right hand side panels. 
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Hubble parameter, Ho, in units of 100 km sec -1 Mpc" 1 , r is the optical depth to reionization 
and n s is the spectral index of the primordial curvature perturbation. As expected from the 
above argument, the effect of Am p is very similar to that of Qb- Due to the larger value of 
Am p , the height of the odd peaks increases compared with that of the even peaks, which is 
a very well-known result for fig 58]. 

Now we discuss how the variations of the parameters change the shape of the CMB power 
spectrum quantitatively. We use four quantities which have been proposed to characterize 
the temperature spectrum well [53]: the position of the first peak £i, the height of the first 
peak relative to the large angular-scale amplitude evaluated at £ = 10, 



AT,, 



2 



the ratio of the second peak £ 2 height to the first, 



AT, 

the ratio of the third peak £3 height to the first, 

*-©)"■ 

where (AT^) 2 = £(£ + l)Cj T /2n. We calculate the response of these four quantities when we 
vary the cosmological parameters. When we vary one parameter, the other parameters are 
fixed. The parameter set we consider is the standard flat-ACDM cosmological parameters 
plus varying physical constants: QBh 2 ,^l m h 2 ,h,T,n s ,a rec ,m e ^ ec ,m Pjrec , where fl m means 
Qb + &dm ■ The fiducial values of the standard cosmological parameters are given in ([2]), and 
those of physical constants equal to the present values. Let us define u>b = ^s^ 2 , uj m = Q m h 2 
and omit the subscript "rec" of physical constants for simplicity, then we have found after 
extensive numerical calculations that the four characteristic measures of the angular power 
spectrum have the following dependence on the cosmological parameters and the physical 
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constants. 



Auih Au m Ah An, Aa Am e Am n 

OA n m AK 1 i QC fi s Q1 7 l 1 AO e l 1(1 1 P 



A£i = 15.1 - 24.0 - 45.3— + 35.6 - 317 + 142 + 19.1 + (6) 

Ub uj m h n s a m e m p 

Auou AuJm Ah An q At 
A#i = 2.84 - 2.80 - 2.02— + 18.3 - 0.64 

uj h uj m h n s r 

+ 4.42^ + 1.05^ + 0.868^, (7) 

a m e m p 

AH 2 = -0.29^ - 0.020^ + 0.53^ + 0.75^ + 0.20^ - 0.351*™* (8) 

Uf, u m n s a m e m p 

Awi AuJrr, An, Aa Am, Am„ , , 

AH 3 = -0.18 + 0.19 + 0.33 + 0.42 0.008 - 0.246 *, (9) 

LUb uj m n s a m e m p 

where the values at the fiducial parameter values are i\ = 219.9, Hi = 6.864, H 2 = 0.4482 
and H 3 = 0.4342, respectively. These formulae are extension of 32] to incorporate the change 
of m p . The coefficients of the AuJb/oJb and Am p /m p are very similar, which agrees with our 
qualitative discussion. Although two parameters and m p have a strong degeneracy, it 
is not exact and we can expect that this small but finite difference actually enables us to 
discriminate these two parameters and give some limit on Am p /m p . More quantitative 
analysis of the degeneracy is performed in the appendix. 



III. THE DILATON DEPENDENCE OF a, m e AND m p 



In this section, we describe our setups for the dilaton dependence of the three phy sical 



constants and relation among them. This generally follows what is adopted in Refs. 
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40] 



As a concrete example, we consider the tree level low energy action of the heterotic string 



in the Einstein frame 
S = 



40 



60 



61 



62 



63 



65[. The action is 



d Xa 



2k 2 m 2 m 



(10) 



where k 2 = 8ttG = %ir/M 2 ^ where M p \ is the Planck scale, $ is the dilaton field, 



is 



an arbitary scalar field, ip is an arbitary fermion, is the gauge covariant derivative 
corresponding to gauge fields with field strength F^ v , and is the gauge field with gauge 
group including SU(3) x SU(2) x U(l). Q is the conformal factor used to transfer from the 
string frame and defined as 



n 



-K$/\/2 



(11) 



Comparing the gauge field strength term — (a'/16K 2 )Q 2 F /lu F fll/ of the action (jTUl) with 
the definition of the Lagrangian density — (l/4g 2 )F /lu F' J ' u where g is the unified coupling 
constant, we get 

1 a'e~^ aS , . 

(12) 



g 2 (M pl ) Ak 2 4« ' 

The gauge coupling constants at low energy scale are calculated using renormalization group 
equations. As for the fine structure constant a, we can check that a does not run practically, 
so a at low energy (or CMB energy scale which is equal to the energy scale at which 
recombination took place) is determined in good approximation as 

a rec « a(M pl ) = = *U^. (13) 

From this expression, our variable Act /a becomes 



Aa _ a mc - a _ y^A* 



- 1 (14) 



a «o 
where A$ = $ rcc - $ . 

From the solution of the one-loop renormalization group equation for the SU(3) coupling 
constant g^ 3 , 

« a _ 9s(M pl ) 2 _ 127T 

4vr 271og(M p VA 2 QCD )- 1 J 

Using this relation and the fact that g(M p i) 2 = 0.1 and that the proton mass m p is propor- 
tional to the QCD energy scale Aqcd, Am p /m p is also expressed through Aa/a as 

Amp ^ m p , rec - m p , now = AA QCD = / 80tt 2 \ _ ^ 

m p m p , now A QCD P ^ 9 1 + As J 

This relation is peculiar in that a small increase in Aa/a results in a large increase in 
Am p /m p , which has a great impact on the parameter estimation. 

Finally, we can read the electron mass term from the action as e K ^^m^^) . Assuming 
that the Higgs vacuum expectation value (H) = const., the time variation of m e is expressed 
as 

Am P m» r «, — m»n „a*/./o /. Aa 



m e m e n V a 



The integration constant is determined by 33 (Aqcd) = 00 
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FIG. 2: One-dimensional marginalized posterior distributions for four models: varying all three 
constants a,m e and m p assuming the relation (fT6|) and (JTTJ) (solid blue curves), varying only two 
constants a and m e assuming the relation (|17p (dot-dashed green curves), varying only one constant 
a (dotted red curves) and varying only one constant m p (dashed magenta curves). 

IV. CONSTRAINTS ON THE VARIATION OF THE COUPLING CONSTANTS 



Based on the model described in the previous section, we constrain the variation of a 
using three kinds of CMB anisotropy spectra, Cf T , Cf , and Cj E of the five-year WMAP 



data 



56 



57|, 



661 ] . Theoretical anisotropy spectra are calculated by the CAMB code [55], |67| 



modified to include the varying physical constants as in Sec. [TTJ We performed the parameter 



estimation using 
CosmoMC code 



Vfarkov-Chain Monte-Carlo (M CMC) techniques implemented in the public 



68 



69]- 



We have run the CosmoMC code on eight Markov chains in each case. To check the 
convergence, we used the "variance of chain means" / "mean of chain variances" R statis- 
tic and adopted the condition R — 1 < 0.03. In all the analysis below, we have also 
incorporated the result of Hubble Key Project of the Hubble Space Telescope (HST) 
on the Hubble parameter Hq, which means that we have imposed Gaussian prior of 
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TABLE I: MCMC results on the mean values and 68% confidence intervals of cosmological pa- 
rameters for four versions: varying a only, varying a and m e , varying a,m e and m p , varying m p 
only. 





varying a only 


varying a and m e 


varying a, m e and m p 


varying m p only 


100tt B h 2 


9 97 +0.10 
z - z ' -0.10 


9 9/1 +0.09 
z - z ^-0.09 


O 1 O + 0.85 

°- lo -0.83 


O 1 q + 0.97 

J - io -0.95 


^DMh 2 


0.109t°;°°J 


o.io8l8:88 6 7 


n ino+0.008 
U.1UO_ 


0.10418:88? 


T 


0-0877±8SS 


0.0864+8:88^ 


n 09 1 6+ 007 


0-092818Z 


n s 


U.JOD_q Q15 


Q61+ - 014 


n 04^+0.020 

U.J4O_ 020 


n 047+0.022 
U - J ^' -0.022 


log(10 10 ^) 


o n6+ - 05 


o n6+ 004 


nn+0.06 
°- uu -0.06 


O m +0.06 

°- ul -0.07 


H 


71.9±?;S 


71.0181 


68.21^ 


69.9l^j 


Aa/a 


-o.ooo89ig;8^ 


-o.ooi9i8:8^ 


U.UU,31_o.0029 




Am p / m p 






-0 22+ - 19 
u - zz -0.19 


-0 19+ - 22 

u - iy -0.23 



Hq = 72 ± 8 km sec -1 Mpc -1 [70|]. In our analysis, we derive the confidence limits for 
standard parameters (Qsh 2 , QoMh 2 , Ho,n s , A s ,t) plus some varying physical constants. 

Figures [2] show ID-marginalized posterior distributions of the eight parameters for the 
case not only a but also m e and m p are changed simultaneously according to the model 
described in the previous section. This is the main result of this paper and for comparison 
with the literature we have also plotted the distributions for two other cases: the case of 
varying only a 35] and the case of varying both a and m e assuming the relation (TlTl) which 
has already been analyzed in 49(. In TABLE U, we also present the mean values and the 68% 
confidence intervals of the cosmological parameters for these three cases. We can conclude 
from these results that the inclusion of the variation of m p results in a tighter limit on the 
variation of a, which is now 

Aa 



8.28 x 1(T 3 < < 1.81 x 1(T 3 , 

a 



(18) 



compared with the limit for the two other cases: —0.024 < Aa/a < 0.019 (for varying 
a, m e ) and —0.028 < Aa/a < 0.026 (for vayring only a), all of which are at 95%C.L. 

In FIG. |2]and TABLE HJ we can also see that, as a cost of tightening of the constraint 
on a, the limit on fish 2 , which is one of the most tightly constrained parameter from CMB 



11 



<3 -0.01 
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0.022 

pB.vech 



2 / \ 3 



Pcrit,0 



a 



FIG. 3: Two-dimensional marginalized posterior distributions of PB^ech 2 / p cr \t,o and Act/a for two 
models: varying all three constants a, m e and m p assuming the relation (|16l) and (|17() (solid blue 
contours) and varying only one constant a (dot-dashed red contours). The plotted pB,rech 2 / p C v\tfl 
(horizontal axis) is actually (pB,rec^ 2 /Pcrit,o)(arec/ao) 3 ; which corresponds to (/2B,o^ 2 //°crit,o)(l + 
Am p /m p ) for the former case and PB,oh 2 / p C rit,o f° r the latter case. In each pair of contours, the 
inner one shows 68%C.L. and the outer one shows 95%C.L. 

data in ordinary analysis, become very loose because of the degeneracy between Qb and 
Am p /m p . This result is reasonable from the argument in Sec. [Til On the other hand, 
since CMB is sensitive to the values of cosmological parameters at the recombination epoch, 
the limit on Qb,?cc m the present scenario is expected to be of the same order as that on 
fls,iec i n the case without m p variation. We can confirm it from FIG. El which shows 2D- 
marginalized posterior distributions of Aa/a-p S]rec /i 2 /p crit for two cases, one with both m p 
and m e variation and the other with no variation in m p and m e . 

One of the reasons why we could obtain tighter constraint than the case without varying 
m p is thought to be the specific relation (Tl6]) in our model. To ensure that statement, we 
have also implemented MCMC calculations in the case that only m p is varied and the other 
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-0.6 -0.4 -0.2 0.2 0.4 -15 -10 -5 



Afi/fj, A<I> x1 ° 3 

FIG. 4: One-dimensional marginalized posterior distributions for two parameters: A////X and A<J> 

two physical constants remain constant, which means that we do not impose any model in 
this case. The result is also shown in FIG. |2] as ID- marginalized posterior distributions. 
Two marginalized curves in the panel showing Am p /m p are very similar and those limits at 
95%C.L. are —0.54 < Am p /m p < 0.36 (for varying m p only) and —0.51 < Am p /m p < 0.17 
(for varying a,m e and m p ), which implies that the dominant factor in constraining the 
variation of physical constants under the relation ( TlTl) and ( fT6j) is Am p /m p and it is the 
specific relation between Aa/a and Am p /m p that determines the order of magnitude of the 
resultant Aa/a limit in the case these parameters are varied simultaneously through the 
dilaton's motion. 

Finally, in our model, we can calculate the limit on not only Aa/a but also A$ and 
A/i/ /i where [i = m p /m e . The results are shown in FIG. HI A$ is the dilaton field variation 



and using the relation ([141) with k = V8nG = (2.43 x 10 18 [GeV])-\ we obtain at 95%C.L. 

- 1.4 x 10 16 [GeV] < A$ < 3.1 x 10 15 [GeV]. (19) 
On the other hand, the constraint on Afi/fi is, at 95%C.L., 

-0.52 < — < 0.17. (20) 
Note that, although it is not so stringent compared with those obtained from other experi- 



ments 37| , constraining this parameter in the recombination epoch is meaningful because in 
the previous works the variation of m p has never been considered. We also expect that future 
CMB measurement such as PLANCK will give more stringent limits of these parameters 
(see Appendix B). 
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V. CONCLUSION 



In summary, we have pointed out that the time variation of Aqcd can affect the CMB 
through the variation in the proton mass. We have shown that there is a notable effect when 
we impose the widely adopted assumptions of unification models found in the literature 
(4^. From the WMAP 5-year data, by MCMC analysis using the CosmoMC code, we have 
obtained the new and stringent constraint on the time variation of the fine structure constant 
as —8.28 x 10~ 3 < Aa/a < 1.81 x 1CT 3 . We have also verified that, in getting the above 
result, the newly considered effect of Am p /m p is very important in that the degeneracy 
between Am p /m p and £lg is tremendously strong and that the small variation of Aa/a 
leads to the large variation of Am p /m p . Furthermore, by including the varying m p , we have 
obtained a limit on the variation of proton-electron mass ratio fi as —0.52 < Afi/fi < 0.17. 

CMB observation will be much more precise to give a lot of information around the 
damping tail of the power spectrum in near future. The difference in the CMB power 
spectra between Am p /m p and Qb appears in the small scale diffusion damping because the 
damping scale is determined by the baryon number density rig, probing the time variation of 
coupling constants including Aqcd by CMB will be even more promising. This is discussed 
more quantitatively in appendix A, and supported by a forecast for future CMB survey in 
appendix B. 
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Appendix A 

In this appendix, we discuss the degeneracies between cosmological parameters using 
the linear algebra. In section [Til we have introduced four variables £±, Hi, if 2, H3 and 
calculated how they depend on the cosmological parameters and the physical constants 
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(ujj,, u) m , h, t, n s , a Iec , m e)Tec , m P!iec ), the results of which are shown in ([6])-([9]). These are 
rewritten in the matrix form as follows. 

Au) m 



AH i 



Hi 

AH 2 
H 2 

\ &Ma I 

\ H 3 / 



( 0.069 -0.109 -0.206 0.162 -1.442 0.646 0.087 \ 

0.414 -0.408 -0.294 2.666 -0.093 0.644 0.153 0.127 

-0.647 -0.045 1.183 1.673 0.446 -0.783 

\ -0.415 0.438 0.760 0.967 -0.018 -0.567 / 



Ah 
h 

An a 
n s 

At 

T 

Aa 

a 

Am e 
m e 
Am, 



V Trip J 

(21) 

where we have divided the both sides of each equation by the fiducial values t\ = 219.9, H% = 
6.864, H 2 = 0.4482, H 3 = 0.4342, respectively. The 4x8 matrix connects the cosmological 
parameters and the CMB variables. Such a rectangular n x m matrix, which we call A, can 
be decomposed in general as 



A = UEV T {U : n x n, X : n x m, V : m x m), 



(22) 



which is known as the singular value decomposition. If n < m and rankA = n, then the 
form of £ is like 



/ 



\ 



Or. 



\ 



I 



(23) 



where <jj is the non-zero singular value of the matrix A and the matrices U and V are 
orthonormal ones. Applying the singular value decomposition to the 4x8 matrix in (I2T]) 
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and using the notation D = [/£, we obtain 

/ -0.672 -1.227 -0.794 0.033 o\ 

2.463 -1.350 0.322 0.008 

2.1145 0.815 -0.488 -0.238 

\ 1.268 0.614 -0.232 0.398 / 

-0.073 -0.037 -0.046 0.788 



D 



( 



V 1 



-0.326 0.209 0.149 

0.478 -0.121 0.067 

-0.028 0.825 -0.042 

0.019 -0.013 0.027 

0.702 0.331 -0.318 

0.0670 0.285 0.817 



-0.024 



-0.018 


0.581 


0.069 


-0.168 


0.029 


0.654 


-0.148 


-0.289 


-0.029 


0.304 


-0.660 


0.473 


-0.003 


-0.253 


-0.419 


-0.166 


0.998 


-0.006 


-0.020 


0.013 


0.018 


0.230 


0.450 


-0.106 


-0.019 


0.049 


0.347 


0.344 


0.019 


0.153 


0.196 


0.712 



Now, OH becomes 



AH i 



Hi 
AH 2 

\ H 3 / 



DV i 



( A^k \ 

Auj m 

Ah 

h 

An a 
n s 
At 

T 

Am v • 
\ m v / 



D 



( 0.788^ + 0.581^ 

n s a 

-0.542^ + 0.654^ 

n s a 

0.478^ - 0.660^ + 0.473^ 

Ufr m e m p 

0.825^ - 0.419^ 
0.999^ 

T 

-0.702^ + 0.451^ 
0.817^ 

h 



(24) 



(25) 



(26) 



\ -0.402^ - 0.447^ + 0.712^ / 

At the last transformation, we have neglected the terms whose coefficients are small and 
irrelevant. Remembering the explict form of the matrix D that the last four columns are 
zero, the last four components in the column vector represent the four degeneracy direc- 
tions, 0.999Ar/r, -0.702Ac; 6 /cj 6 + 0.451 Am e /m e , 0.817 Ah/h, -0A02Au b /u b -0.U7Ah/h+ 
0.712Am p /m p . 

In the above case, where all the variations of physical constants are included in the 
analysis, the final results are complicated and it is difficult to understand the degeneracies 
cleary. Therefore, we consider simpler cases in order now. 
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First, we treat only standard cosmological parameters, so (|2"T|) is reduced to 



11 

AH i 
AH 2 

\ Ma I 

\ h 3 / 



( 0.069 -0.109 -0.206 0.162 \ 

0.414 -0.408 -0.294 2.666 -0.093 

-0.647 -0.045 1.183 

\ -0.415 0.438 0.760 



Aw m 

Mm 

Ah 
h 

n a 

V ¥ / 



(27) 



The singular value decomposition of the 4x5 matrix gives 



D 



( -0.190 0.131 -0.052 -0.172 o\ 

-2.713 0.420 0.039 0.017 

-1.171 -0.626 -0.239 0.006 

\ -0.708 -0.609 0.261 -0.031 / 



and we get the approxmated relation 

/ 4^ \ 



V 



CI 

AH i 

H! 

AH 2 
H 2 

\ AH3 I 

\ H 3 / 



DV l 



( \ 




Au> m 








Ah 
h 


« D 


An s 








\ ¥ J 





( 



-0.009 


0.094 


0.090 


0.879 - 


-0.443 


-0.157 


0.452 


0.888 


-0.005 


0.141 - 


-0.075 


0.982 


0.056 





0.043 


-0.991 


An 3 
n s 


\ 



-0.991 0.027 \ 

-0.066 -0.041 

0.0787 -0.028 

0.080 -0.053 

0.031 0.997 ) 

(28) 



0.879^ - 0.443^ 
0.452^ + 0.888^ 



0.983^ 

h 



\ 



0.997 



At 



(29) 



/ 



This tells us that the four standard cosmological parameters, n s , u^, u; m and h can be deter- 
mined by observing the four variables £i,Hi, H 2 and if 3, however the optical depth r cannot 
be limited sufficiently, which is the direct consequence of the fact that we cannot tightly 
constrain r only through the temperature anisotropy spectrum. 

In the next case, we add one physical constant Aa/a to the parameter vector; 

/ A£i> \ 

/ 0.069 -0.109 -0.206 0.162 -1.442 \ 
0.414 -0.408 -0.294 2.666 -0.093 0.644 
-0.647 -0.045 1.183 1.673 
\ -0.415 0.438 0.760 0.967 / 



AH 1 

AH 2 
H 2 

\ Ml I 

\ H 3 / 



Au m 

Ah 
h 

An s 
n s 

At 



(30) 
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From the singular value decomposition and the relevant approximation, we obtain 

/ -0.694 -1.219 -0.436 -0.094 o\ 
-1.2123 0.172 0.045 



D 



2.540 
1.970 
\ 1.182 

0.062 



0.810 
0.541 

-0.044 

-0.342 0.212 0.156 

0.789 -0.459 0.106 
0.500 0.861 
0.056 
\ 0.062 0.029 



-0.162 -0.239 
-0.354 0.246 / 

-0.049 0.810 -0.019 0.579 \ 

-0.528 0.029 0.731 

-0.200 -0.043 0.337 

0.069 -0.032 0.025 

0.032 0.997 -0.001 

-0.141 0.052 0.125 ) 



(31) 



0.048 
0.051 
-0.978 



(32) 



and 



A.Hi 



AH 2 
H 2 

\ ASA I 

\ H 3 / 



DV 1 



( A^k \ 




Auj m 








Ah 




h 
An 3 


« D 


n a 




At 




T 




V ^ > 





0.810^ + 0.579^ 

n a a 

-0.342^ - 0.528^ + 0.731^ 
> A< ^ - 0.459^ + 0.337^ 



0.789- 



0J b 



\ 



0.500^ + 0.861^ 
0.997^ 



-0.978^ 



(33) 
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The last equation ( |33l) makes a contrast with ( !29l . The newly included parameter Aa/a 
is now determined together with n s , uuj, and u m , on the other hand the Hubble parameter 
h becomes the unlimited parameter like r. This is exactly what is expected because in the 



previous paper [35|J we have made sure that, in constraining Aa/a, the HST prior on the 
Hubble parameter is very useful, which implied that including Aa/a in the analysis weakens 
the constraint on h. 

As the third case, which we are most interested in, we add the new parameter Am p /m p ] 

( A£b \ 

( 0.069 -0.109 -0.206 0.162 0.087 \ 
0.414 -0.408 -0.294 2.666 -0.093 0.127 
-0.647 -0.045 1.183 -0.783 
\ -0.415 0.438 0.760 -0.567 / 



( ^\ 

AH i 

i?i 
AH 2 

H 2 

\ AEa I 

\ H 3 / 



Au m 

Ah 
h 

An s 
n s 

At 



Am 



(34) 



\ m p / 
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Similar decomposition and approximation lead to 

/ -0.177 0.170 -0.056 0.172 o\ 
-2.676 0.629 0.0305 



D = 



0.017 

-1.260 -0.890 -0.230 -0.006 

\ -0.776 -0.763 0.281 0.031 0/ 
/ 



(35) 



V s 



0.002 


0.088 


0.088 


-0.985 


0.027 


0.114 


0.647 


-0.316 


-0.122 


0.040 


-0.033 


0.681 


0.301 


0.936 


0.0190 


0.102 


-0.0210 


0.145 


-0.139 


0.075 


-0.983 


-0.080 


0.053 


-0.002 


0.0380 


0.003 


0.046 


0.033 


0.997 


0.019 


-0.685 


0.098 


0.095 


0.098 


0.004 


0.709 



(36) 



and 



AHi 
Hi 

AH 2 



H 2 



DV 1 



/ 




\ 














Au m 
















Ah 








h 




« D 




An s 








n s 








At 








T 






V 


Am p 


J 




m p 





( -0.985^ \ 

-0.647^ + 0.681^ 

ui b m p 

0.936^ 
-0.983^ 

h 



0.997^ 



\ -0.685^ 



(37) 



UJf, 



0.709^ / 

m v / 



It is found that there are two degeneracy directions and one of them is — 0.685Ac^/c<;fc + 
0.709Am p /m p . This means that Ub and m, p are strongly degenerate, which could be naively 
expected from just comparing the coefficients of the two parameters Acub/cui, and Am p /m p . 

Now, we include the small scale information in the analysis to study how the damping 
tail is crucial for parameter determination from CMB in the third case. We make use of Tf 4 
and if 5 , which is defined in the same way as H 2 ,H 3 , that is, 



Ha = 



AT,, 
AT, 



Hz 



ATV 



AT/ 



(38) 



In other words, H 4 is the ratio of the fourth peak £ 4 to the first peak i x height and H 5 is the 
ratio of the fifth peak £ 5 to the first peak £±. Adding the relations between ATf 4 / T/ 4 , AH 5 /H 5 
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and the variations of the cosmological parameters, (1341) becomes 



/ \ 

AH i 



AH 2 

AH 3 

H 3 
AH 4 

H 4 

\ Ms I 

\ H s / 



( 0.069 -0.109 -0.206 0.162 

0.414 -0.408 -0.294 2.666 

-0.647 -0.045 1.183 

-0.415 0.438 0.760 

-0.478 0.315 1.608 

\ 0.665 0.493 1.828 



0.087 
-0.093 0.127 
-0.783 
-0.567 
-0.868 
-0.794 





/ A<^ \ 
















Ah 




h 




An s 




n s 




At 

T 




Am p 
\ m p / 



(39) 



We can recast this equation through the relevant approxmation in 



/ Ah \ 

h 
AH! 



Hi 

AH 2 
H 2 

AHa 
H 3 

AH 4 
H 4 



( 0.1378 -0.198 0.085 -0.012 -0.165 0.001 \ 

2.528 -1.024 0.343 -0.042 0.023 -0.001 

1.332 0.680 0.409 0.170 -0.013 -0.003 

0.897 0.648 -0.027 -0.194 -0.014 -0.004 

1.785 0.684 0.104 -0.0469 0.012 0.006 

\ 2.007 -0.046 -0.789 0.069 -0.013 -0.001 / 



0.960 



An. 



-0.728^ 





0.658 Amp 






Aujt 


0.579 AuJm 


iv b 




Au) m 


0.582 Amp 







0.959^ 



V 



0.994^ 



(40) 

All the elements of the last column in 6 x 6 matrix, which correspond to the coefficients 
of the parameter r, are nearly zero, which is almost the same situation occurred in the 
previous case without Aif 4 /iJ 4 and AH 5 /H 5 . Unlike the previous case, there are no other 
(nearly) zero columns, which states that, as mentioned in section |V] the strong degeneracy 
between uib and m p can be resolved by observing not only £i, Hi, H 2 , H 3 but also H^, H 5 
which possess the damping tail information. 



Appendix B 

In this appendix, we briefly report a forecast for constraints from the future CMB sur- 



vey, in particular, the PLANCK survey 71|. The parameters for the instrumental design 
for PLANCK survey are shown in TABLE [III arid mock data are generated by the pub- 
licly available FuturCMB code [3]. The fiducial parameter values are the same as those 
adopted in section [Til and we perform the MCMC analysis in the same way as in section 
IIVI Throughout this appendix, we assume the model between time variations of physical 
constants introduced in section IIII1 
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TABLE II: Assumed instrumental parameters for PLANCK survey. f s k y is the fraction of the ob- 
served sky, #fwhm is Gaussian beam width at FWHM, ut and up are temperature and polarization 
noise, respectively. 



fsky 


bands [GHz] 


$fwhm [arcmin] 




ap[fjK] 


0.65 


100 


9.5 


6.8 


10.9 




143 


7.1 


6.0 


11.4 




217 


5.0 


13.1 


26.7 



TABLE III: MCMC results on the mean values, 68% and 95% confidence intervals of cosmological 
parameters. 





varying a, m e and m p 


won B h 2 


9 97/1+0.025+0.050 
z - z '*-0.025-0.049 


tt DM h 2 


n i i i +0.002+0.003 
u - iii -0.002-0.003 


T 


n ns7+ - 002 + - 009 

U.UO^_ 003-0 .009 


n s 


n qrq+0.009+0.017 
U - JUO -0.009-0.018 


log(10 10 A fl ) 


q nR c+0.011+0.023 
°- uod -0.011-0.023 


H 


71 q+1. 0+2.1 
a - y -l. 1-2.0 


Aa/a 


o nnooo2i +0 - 0002+0 - 0003 

U.UUUUUZJ._ .0002-0.0003 


Am/a* 


n nnnQ+°- 016 + - 030 
u - uuuo -0. 015-0.031 



In TABLE IHIl we show the forecast for the constraints on the cosmological parameters 
as the mean values, 68% and 95% confidence intervals. As is commented in the concluding 
section and appendix A, future CMB survey can constrain time variation of the physical 
constants more tightly and the limit on Aa/a or A///// from PLANCK experiment is ex- 
pected to be about one order tighter than the constraint from the WMAP 5-year data and 
comparable to the constraint from BBN 
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